Electroproduction of Charmonia 
off Protons and Nuclei 



Yu.P. Ivanov^, B.Z. Kopeliovich***, A.V. Tarasov f and J. Hiifner* 



(N 
O 
O 



* Institut filr Theoretische Physikder Universitat, 69120 Heidelberg, Germany 
' Joint Institute for Nuclear Research, Dubna, 141980 Moscow Region, Russia 
| ** Institut filr Theoretische Physik der Universitat, 93040 Regensburg, Germany 

^Max-Planck Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 

<D 

£^ ' Abstract. Elastic virtual photoproduction of charmonia on nucleons is calculated in a parameter 

free way with the light-cone dipole formalism and the same input: factorization in impact parame- 
ters, light-cone wave functions for the photons and the charmonia, and the universal phenomenolog- 
ical dipole cross section which is fitted to other data. The charmonium wave functions are calculated 
, with four known realistic potentials, and two models for the dipole cross section are tested. Very 

good agreement with data for the cross section of charmonium electroproduction is found in a wide 
£S| ! range of s and Q 2 . Using the ingredients from those calculations we calculate also exclusive elec- 

ts) ■ troproduction of charmonia off nuclei. Here new effects become important, (i) color filtering of the 

cc pair on its trajectory through nuclear matter, (ii) dependence on the finite lifetime of the cc fluc- 
tuation (coherence length) and (iii) gluon shadowing in a nucleus compared to the one in a nucleon. 
Total coherent and incoherent cross sections for C, Cu and Pb as functions of s are presented. The re- 
sults can be tested with future electron-nucleus colliders or in the peripheral collisions of relativistic 
heavy ions. 
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C^: INTRODUCTION 
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In contrast to hadro-production of charmonia, where the mechanism is still debated, 
electro(photo)production of charmonia seems better understood: the cc fluctuation of 
the incoming real or virtual photon interacts with the target (proton or nucleus) via the 
dipole cross section o qq and the result is projected on the wave function of the observed 
hadron [1]. The aim of this paper is not to propose a conceptually new scheme, but to 
calculate within this approach as accurately as possible and without any free parameters. 
Wherever there is room for arbitrariness, like form for the color dipole cross section and 
for for charmonium wave function, we use and compare other author's proposals, which 
have been tested on different data. 

In the light-cone (LC) dipole approach the virtual photoproduction of charmonia (here 
*F stands for J/y or looks as shown in Fig. 1 [1]. The corresponding expression for 
the forward amplitude reads 

l 

s,Q 2 ) = Lj da J d%^\a,r T )a qq (r T: s)^(a,r T ,Q 2 ). (1) 
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Here the summation runs over spin indexes /i, jli of the c and c quarks, Q 2 is the photon 



a qq (r T ,s) 




FIGURE 1. Schematic representation of the amplitude for the reaction yp — > *¥p in the rest frame of 
the proton. The cc fluctuation of the photon with transverse separation r T and cm. energy y/s interacts 
with the target proton via the dipole cross section O q Jr T ,s) and produces a x ¥. 



virtuality, ^(oc, r T , Q 2 ) is the LC distribution function of the photon for a cc fluctuation 
of separation r T and relative fraction a of the photon LC momentum carried by c or 
c. Correspondingly, ^(a,?^) is the LC wave function of J/xj/ or y/. The dipole cross 
section <J qq -(r T ,s) mediates the transition. 

This paper is organized as follows. First we review the expressions of the factorized 
LC approach to electroproduction of heavy quarkonia (more details one can find in [2]) 
and then compare with available experimental data for J / y production. The calculations 
which are parameter free demonstrate good agreement with data. 

After it we analyse exclusive electroproduction of charmonia off nuclei yA — > ^X, 
where X = A (coherent) or X = A* (incoherent, where A* is an exited state of the 
A-nucleon system). In these processes new phenomena are to be expected, (i) Color 
filtering, i.e. inelastic interactions of the cc pair on its way through the nucleus is 
expected to lead to a suppression of *P production relative to Ao ^, p . (ii) Production 
of a cc pair in a nucleus and its absorption are also determined by the values of the 
coherence length (lifetime of the cc fluctuation), 

/ - 2V r, 2V (2) 
lc Q 2 +M?, Q 2 +M 2 /v {l) 

where v is the energy of the virtual photon in the rest frame of the nucleus, (iii) Since 
the dipole cross section o qq - also depends on the gluon distribution in the target (p of 
A), nuclear shadowing of the gluon distribution is expected to reduce o qq in a nuclear 
reaction relative to the one on the proton. 

The predictions in this paper may be used for planning of future experiments for 
electron-nucleus collisions at high energies like in the eRHIC project. Another possi- 
bility to observe photoproduction off nuclei is heavy ion relativistic collisions (see, for 
example, review [4]). In the last section we present our results for J/\\f production in 
such processes. 



LIGHT-CONE DIPOLE APPROACH 

The LC variable describing longitudinal motion of the quarks is the fraction a = / Py 
of the photon LC momentum p+ = E y + p y carried by the quark or antiquark. a 



is Lorentz-boost invariant. In the nonrelativistic approximation (assuming no relative 
motion of c and c) a = 1/2 (e.g. [1]), otherwise one should integrate over a (see Eq. (1)). 



Photon wave function 

For transversely (T) and longitudinally (L) polarized photons the perturbative photon- 
quark distribution function in Eq. (1) reads [5, 6], 

^f(a^Q 2 ) = ^^Z cX ^d T ^K (er T ) , (3) 

where % d = io y x^, X an d X are the spinors of the c-quark and antiquark respectively; 
Z c = 2/3. K (er T ) is the modified Bessel function with e 2 = a(l - a)Q 2 + m 2 .. The 
operators TL have the form: 

T = m c a-e r +i(l-2a)(a-n)(e r -V rT ) + (nxe r )-V rT , (4) 
d L = 2Qa(l-a)o-n , (5) 

where n = p/p is a unit vector parallel to the photon momentum and e is the polarization 
vector of the photon. Effects of the non-perturbative interaction within the qq fluctuation 
are negligible for the heavy charmed quarks. 



Charmonium wave function 

The charmonium wave function is well defined in its rest frame where one can rely on 
the Schrodinger equation. As soon as the rest frame wave function is known, one may 
be tempted to apply the Lorentz transformation to the cc pair as it would be a classical 
system and boost it to the infinite momentum frame. However, quantum effects are 
important and in the infinite momentum frame a series of different Fock states emerges 
from the Lorentz boost. Therefore the lowest \cc) component in the infinite momentum 
frame does not represent the \cc) in the rest frame. We rely on the widely used procedure 
[1 1] for the generation of the LC wave functions of charmonia. 

In the rest frame the spatial part of the cc pair wave function satisfying the Schrodinger 
equation 

(-^ + V^V nl Jr)=E nl V nlm (r) (6) 
is represented in the form 

nr)=^ nl {r)-Y lm {e^), (7) 



where r is 3-dimensional cc separation, ^^(r) and Y, (6, cp) are the radial and orbital 
parts of the wave function. The following four potentials V (r) have been used 



'COR": Cornell potential [12], 



V (,■) = — + 4 (8) 



r 

with = 0.52, a = 2.34GCV" 1 and m c = 1.84GeV. 

• "BT": Potential suggested by Buchmiiller and Tye [13] withm c = 1.48 GeV. It has a 
similar structure as the Cornell potential: linear string potential at large separations 
and Coulomb shape at short distances with some refinements, however. 

• "LOG": Logarithmic potential [14] 

V(r) = -0.6635 GeV + (0.733 GeV) log(r • 1 GeV) (9) 

with m c = 1.5 GeV. 
. "POW": Power-law potential [15] 

V(r) = -8.064GeV + (6.898 GeV) (r-1 GeV) 01 (10) 

withm c = 1.8 GeV. 

The shapes of the four potentials differ from each other only at large r (> 1 fm) and very 
small r (< 0.05 fm) separations. Note, however, that COR and POW use m c xs 1.8 GeV, 
while BT and LOG use m c ^1.5 GeV for the mass of the charmed quark. This difference 
will have significant consequences. 

For the ground state 15 all the potentials provide a very similar behavior for the radial 
part ^* n /(r) at r > 0.3 fm, while for small r the predictions differ by up to 30%. The 
peculiar property of the 25 state wave function is the node at r m 0.4 fm which causes 
strong cancellations in the matrix elements Eq. (1) and as a result, a suppression of 
photoproduction of y/ relative to J/\j/ [1, 16]. 

The lowest Fock component \cc) in the infinite momentum frame is not related by 
simple Lorentz boost to the wave function of charmonium in the rest frame. This makes 
the problem of building the LC wave function for the lowest \cc) component difficult, no 
unambiguous solution is yet known. There are only recipes in the literature, a simple one 
widely used [1 1], is the following. One applies a Fourier transformation from coordinate 
to momentum space to the known spatial part of the non-relativistic wave function (7), 
m(r) =^ ^(P), which can be written as a function of the effective mass of the cc, 
M 2 = 4(p 2 + m 2 ), expressed in terms of LC variables 

In order to change integration variable p L to the LC variable a one relates them via 
M, namely p L = (a — \/2)M{p T ,a). In this way the cc wave function acquires a 
kinematical factor 

/ 2, 2\3/4 

np) * ^ (p 2 +m 2 y/ 2 ' X¥{a,pT) " * ¥iaM ' (12) 



This procedure was used in [17] and the result is applied to calculation of the ampli- 
tudes (1). The result was discouraging, since the y/ to J/\j/ ratio of the electroproduction 
cross sections are far too low in comparison with data. However, the oversimplified 
dipole cross section o qq -(r T ) oc r \ was used, and what is even more essential, the impor- 
tant ingredient of Lorentz transformations, the Melosh spin rotation, was left out. 

The 2-dimensional spinors % c and % s describing c and c respectively in the infinite 
momentum frame are known to be related via the Melosh rotation [1 1, 18] to the spinors 
Xc and %c m me rest frame: ~%~ c = R(a,p T ) Xc and = R(l - a, -p T ) x- v where the 
matrix R(a,p T ) has the form: 

R(a,p T ) = ^ + aM -noxn]p T ^ 
'(m c + aM) 2 + p 2 



Since the potentials we use contain no spin-orbit term, the cc pair is in 5-wave. In this 
case spatial and spin dependences in the wave function factorize and we arrive at the 
following LC wave function of the cc in the infinite momentum frame 

®^Xa,p T ) = U^Xa,p T )-<S> v (a,p T ) , (14) 

where 

U^\a,p T ) = x^R\a,p T ) a ■ e v o y R*(l - a, -p T ) a y 7 l % ■ (15) 

Now we can determine the LC wave function in the mixed longitudinal momentum - 
transverse coordinate representation: 

^\a,r T ) = ^J d 2 p T e- l ^r^)(aJ T ) . (16) 



Phenomenological dipole cross section 

The color dipole cross section o q q(r T ,s) is poorly known from first principles. It is 

expected to vanish <x r\ lnr r at small r T ^0 due to color screening [7] and to level off at 
large separations. We use a phenomenological form which interpolates between the two 
limiting cases of small and large separations. Few parameterizations are available in the 
literature, we choose two of them which are simple, but quite successful in describing 
data and denote them by the initials of the authors as "GBW" [8] and "KST" [9]. 



"GBW": oJr T ,x) = 23.03 



mb , (17) 



1 _ e -4AoW 

/ n0.144 

r Q (x) = 0.4 fm, 

where x Q = 3.04 ■ 10~ 4 . The proton structure function calculated with this parameteriza- 
tion fits well all available data at small x and in wide range of Q 2 [8]. However, it obvi- 
ously fails describing the hadronic total cross sections, since it never exceeds the value 



23.03 mb. The ^-dependence guarantees Bjorken scaling for DIS at high Q 2 , however, 
Bjorken x is not a well defined quantity in the soft limit. Instead we use the prescription 
of [10], x = (M 2 + Q 2 )/s, where My is the charmonium mass. 

This problem with limited dipole cross section as well as the difficulty with the 
definition of x have been fixed in [9]. The dipole cross section is treated as a function of 
the cm. energy \fs, rather than x, since \fs is more appropriate for hadronic processes. 
A similarly simple form for the dipole cross section is used 



KST": o qq (r T ,s) = o Q (s) 



1 



(18) 



The values and energy dependence of hadronic cross sections are reproduced with the 
following expressions 



23.6 — 



0.08 



3 rl(s) 



mb , 



r (s) = 0.88 — fm . 



-0.14 



(19) 
(20) 



The energy dependent radius r Q (s) is fitted to data for the proton structure function 
FP(x,Q 2 ), s Q = 1000 GeV 2 and the mean square of the pion charge radius (r 2 h ) = 

0.44 fm . The improvement at large separations leads to a somewhat worse description 
of the proton structure function at large Q 2 . Apparently, the cross section dependent on 
energy, rather than x, cannot provide Bjorken scaling. Indeed, parameterization (18) is 
successful only up to Q 2 « 10 GeV 2 . 

In fact, the cases we are interested in, charmonium production and interaction, are just 
in between the regions where either of these parameterization is successful. Therefore, 
we suppose that the difference between predictions using Eq. (17) and (18) is a measure 
of the theoretical uncertainty which fortunately turns out to be rather small. 



ELECTROPRODUCTION OFF PROTONS 

Having Eq. (1) and the expressions from the previous section (LC wave functions and 
dipole cross section), we can calculate cross sections for the virtual photoproduction 

fp — > 

/.^_ l^(^6 2 )l 2 + £|^6 2 )l 2 nu 

where e is the photon polarization (for HI data (e) = 0.99); B is the slope parameter in 

reaction f*P^ VP- We use the experimental value [20] 5 = 4.73 GeV" 2 . ^ T L includes 
also the correction for the real part of the amplitude: 

~ , / it dln^ TT (s, Q 2 )\ 
j# t £ S , Q 2 ) = ^ TX (s, Q 2 )[\-i^ ^ ) , (22) 



s 



where we apply the well known derivative analyticity relation between the real and 
imaginary parts of the forward elastic amplitude [19]. The correction from the real part 
is not small since the cross section of charmonium electroproduction is a rather steep 
function of energy. 
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FIGURE 2. Integrated cross section for elastic photoproduction with real photons (Q 2 = 0) calculated 
with GBW and KST dipole cross sections and for four potentials to generate J/y wave functions. 
Experimental data points from the HI [20], E401 [21], E516 [22] and ZEUS [23] experiments. 
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FIGURE 3. The ratio of to J / \j/ virtual photoproduction cross sections as a function of the photon 
virtuality Q 2 at energy y/s = 90 GeV for four potentials and with GBW and KST parameterizations for 
the dipole cross section. Experimental data points from the HI experiment [24]. 

We present here only results for energy dependence of the J/y/ integrated cross section 
(Fig. 2) and for Q 2 dependence of the ratio y/ to J/\j/ photoproduction (Fig. 3). More 
results are presented in [2]. All calculations are performed with GBW and KST parame- 
terizations for the dipole cross section and for wave functions of the J/\j/ calculated from 
BT, LOG, COR and POW potentials. One can see that there are no major differences for 



the results using the GBW and KST parameterizations. The BT and LOG potentials de- 
scribe the data very well, while the potentials COR and POW underestimate them by 
a factor of two. The different behavior has been traced to the following origin: BT and 
LOG use m c « 1.5 GeV, but COR and POW m c « 1.8 GeV. While the bound state wave 
functions of J/\j/ are little affected by this difference, the photon wave function Eq. (3) 
depends sensitively on m c . 

It turns out that the effects of spin rotation have a gross impact on the ratio R = 
a{Y)/a{J/^f). This effects add 30-40% to the J/\j/ electroproduction cross section. But 
they have a much more dramatic impact on i// increasing the electroproduction cross 
section by a factor 2-3. This spin effects explain the large values of the ratio R observed 
experimentally. Our results for R are about twice as large as evaluated in [25] and even 
more than in [17]. 



ELECTROPRODUCTION OFF NUCLEI 

Exclusive charmonium production off nuclei, yA — > *¥X is called coherent, when the 
nucleus remains intact, i.e. X = A, or incoherent, when X is an excited nuclear state 
which contains nucleons and nuclear fragments but no other hadrons. The cross sections 
depend on the polarization £ of the virtual photon (in all figures below we will imply 

£ = 1), 



<7 yA ( Sl Q z ) = a T (j, <2 ) + £ o L (s, Q z ) , 



(23) 



where the indexes T,L correspond to transversely or longitudinally polarized photons, 
respectively. At high energies the coherence length Eq. (2) may substantially exceed 
the nuclear radius (l c 3> R A ). In this case the transverse size of the cc wave packet is 
"frozen" by Lorentz time dilation, i.e. it does not fluctuate during propagation through 
the nucleus, and the expressions for the incoherent (inc) and coherent (coh) cross sec- 
tions are simple [1]: 



°¥£(s,Q z ) = 



d 2 b 



Jt c T °Hs,Q\b) 



4 n i(^Q 2 ) = Jd 2 b\^ n i(s,Q 2 ,b) 



\6nB(s) ' 



(24) 
(25) 



where T A (b) = J^dzp^b.z) is the nuclear thickness function given by the integral of 
the nuclear density along the trajectory at a given impact parameter b and expressions 

for jZ{s,Q 2 ,b) are given by Eqs. (22) and (1) with replacement 



1 — exp 



-° qq :(r T , S )T A (b)/2 
-a qq (r Tl s)T A (b)/2 



coh. 



inc. 



(26) 



But for charmonium production off nuclei the "frozen" approximation is not enough 
and additional phenomena should be taken into account: effects of finite coherence 
length and gluon shadowing. 



Finite coherence length 



The "frozen" approximation (26) is valid only for l c 3> R A and can be used only at 
y/s > 20 30 GeV. The low-energy part should be corrected for the effects related to 
the finiteness of l c . A strictly quantum-mechanical treatment of a fluctuating qq pair 
propagating through an absorptive medium based on the LC Green function approach 
has been suggested recently in [26] . However, an analytical solution for the LC Green 
function is known only for the simplest form of the dipole cross section o qq (r T ) oc r \. 
With a realistic form of a qq (r T ) it is possible only to solve this problem numerically, 
what is still a challenge. Here we use the approximation suggested in [27] to evaluate 
the corrections arising from the finiteness of l c by multiplying the cross sections for 
coherent and incoherent production evaluated for l c — > °o by a kind of formfactor F coh 
and F inc respectively: 

a A {s,Q 2 ) =s> a A {s,Q 2 ) ■ F (sJ c (s,Q 2 )) , 



where 



■pcoh 



(S,lc) 



d 2 b 



F inc (s,l c ) 
F 2 (s,b,z,l c ) 



dz9 A {b,z)F x {s,b,z)^ /(...)| /c=00 , 
Jd 2 b J dzp A (b,z) \F l (s,b,z)-F 2 (s,b,zJ c )\ 2 /(. 

CO 

( 1 1 , / 

ex P -2°vn( s ) J d zP A ( b ^') 

z 

^ N (s) J dz'p^F^e-^-^. 



(27) 

(28) 
l/ c =oo, (29) 
(30) 
(3D 



For the charmonium nucleon total cross section (J w (s) we use our previous results [2]. 



Gluon shadowing 

The gluon density in nuclei at small Bjorken x is expected to be suppressed compared 
to a free nucleon due to interferences. This phenomenon, called gluon shadowing, 
renormalizes the dipole cross section, 

o qq (r T ,x)^o qq (r T ,x)R G (x,Q 2 ,b). (32) 

where the factor R G (x, Q 2 ,b) is the ratio of the gluon density at x and Q 2 in a nucleon 
of a nucleus to the gluon density in a free nucleon. No data are available so far which 
could provide information about gluon shadowing. Currently it can be evaluated only 
theoretically. To calculate function R G we use the approach developed in [9] and applied 
for charmonia production in [3]. 



Numerical results 



Combining the results above (i.e. including finite coherence length and gluon shad- 
owing) we obtain the results for charmonia (J/y/ and yr') electroproduction off nuclei. 
As it is common practice we express nuclear cross sections in the form of the ratio 



Ac yp {s,Q 2 ) 



(33) 



where the numerator stands for the expression Eq. (23) (with Eq. (24) for coherent 
and Eq. (25) for incoherent scattering) and the denominator is given by Eq. (21). We 
present here only s dependences (plots for Q 2 dependences and momentum transfer k T 
distributions are given in [3]) for coherent (Fig. 4) and incoherent (Fig. 5) production 
of charmonia with the GBW [8] parameterization for the dipole cross section. KST [9] 
parameterization gives quite close results (they differ at most 10% at high energies). 
It is not a surprise that the ratios for coherent production exceed one: in the absence 
of cc attenuation the forward coherent production would be proportional to A 1 , while 
integrated over momentum transfer it behaves as A 4 / 3 . It is a result of our definition 
Eq. (33) that R^ h exceeds one. 
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FIGURE 4. Ratios R"'^ and R c °, h for coherent production on carbon, copper and lead as a function of 

y/s and at Q 2 = calculated with GBW parameterization of a qC] . Solid curves display the modifications 
due to the finite coherence length l c and gluon shadowing corrections while the dashed lines are without 
("frozen" approximation). 

We see that in comparison with the "frozen" approximation, l c corrections noticeable 
change ratios at low energies (especially for i//). For coherent production in the limit 
of low energies l c — > the strongly oscillating exponential phase factor in (28) makes 
the integral very small and thus F coh w 0. Then the cross section rises with l c unless it 
saturates at l c ^> R A when the phase factor becomes constant. Apparently, this transition 
region is shifted to higher energies for larger nuclear radius. For incoherent production 
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FIGURE 5. Ratios R'^' for J/y and y/ 1 incoherent production on nuclei as function of y/s and at Q 2 = 0. 
The meaning of the different lines is the same as in Fig. 4. 



the observed nontrivial energy dependence is easy to interpret. At low energies l c <C R A 
and the photon propagates without any attenuation inside the nucleus where it develops 
for a short time t c a cc fluctuation which momentarily interacts to get on mass shell. The 
produced cc pair attenuates along the path whose length is a half of the nuclear thickness 
on the average. On the other hand, at high energies when l c ^> R A the cc fluctuation 
develops long before its interaction with the nucleus. As a result, it propagates through 
the whole nucleus and the mean path is twice as long as at low energies. This is why the 
nuclear transparency drops when going from the regimes of short to long l c . 

At high energies gluon shadowing becomes important. We see that the onset of gluon 
shadowing happens at a cm. energy of few tens of GeV. Remarkably, the onset of 
shadowing is delayed with rising nuclear radius. Nuclear suppression of production 
becomes stronger with energy. This is an obvious consequence of the energy dependence 
of o q q(r T ,s), which rises with energy. For the suppression is rather similar as for 
J/xy. In particular we do not see any considerable nuclear enhancement of i// which has 
been found earlier [1, 28], where the oversimplified form of the dipole cross section, 
o Jr T ) oc r\ and the oscillator form of the wave function had been used. 



HEAVY ION PERIPHERAL COLLISIONS 

The large charge Z of heavy nuclei gives rise to strong electromagnetic fields: the photon 
field of one nucleus can produce a photo-nuclear interaction in the other. The cross 
section of the charmonia photoproduction by the induced photons reads 

k^-= [d 2 b fd 2 b'n{k,V -b)o A (b,s), (34) 



where k is the photon momentum. Photon flux induced by projectile nucleus with Lorenz 
factor 7 is 

, - N a m Z 2 k 2 7 ( bk\ 
n(k 1 b) = -^ r Kl^-y (35) 



Cross sections o A (b, s) for coherent and incoherent production are 

2 



(s,b) 



?? h (s,b) 



T{b) 
16kB(s) 



(36) 
(37) 



where expressions M{s 1 b) correspond to M in Eqs. (24, 25) at Q 2 = 0. Our predictions 
for coherent J/yr production at RHIC and LHC energies are presented on Fig. 6. We see 
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FIGURE 6. Rapidity distribution for coherent J/y production in heavy ion peripheral gold-gold colli- 
sions at RHIC (left) and lead-lead collisions at LHC (right) calculated with the GBW parameterization of 
<7^. Solid curves display the modifications due to gluon shadowing and finite coherence length l c while 
the dashed lines are without ("frozen" approximation). 

that l c corrections modify the distribution at the edges (positive and negative y) while 
suppression around y = is provided mostly by gluon shadowing (especially for LHC, 
where energies should be much higher). 



CONCLUSION 

In this paper we use the LC dipole approach for description of charmonium electropro- 
duction off protons and nuclei. We have no free parameters and our calculations are in 
good agreement with existing experimental data. Our predictions can be tested in future 
experiments at high energies with electron-nuclear colliders (eRHIC) and in peripheral 
heavy ion collisions (RHIC and LHC). 
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